Transport properties of a spin-| Heisenberg chain with an embedded spin-S* impurity 
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of the impurity by the chain eventually leading to the cutting or healing of the host chain. Numerical 
results are supported by analytical arguments obtained in the strong host-impurity coupling regime. 
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I. INTRODUCTION 

The unconventional thermal transport properties of 
low dimensional quantum magnets have drawn the at- 
tention of the condensed matter society. 1-15 In particu- 
lar, the transport properties of spin chain materials are 
described by the one-dimensional (ID) spin-g Heisen- 
berg model, where the large exchange coupling along 
the chain 16-18 accounts for an extraordinary high and 
anisotropic thermal conductivity. 19 Actually, the pure 
anisotropic Heisenberg model (AHM) was shown to ex- 
hibit ballistic thermal transport at any temperature T, 
attributed to its integrability. 20-22 In accord with the the- 
oretical predictions, evidence of ballistic thermal trans- 
port in real materials has been consolidated lately using 
samples of very high purity. 23 However, defects may have 
a dramatic effect on the transport properties of these sys- 
tems due to the reduced dimensionality. 

In this work, we try to shed light on the effect of a 
single, spin-S* magnetic impurity, on the transport prop- 
erties of a spin-| Heisenberg chain. We address the be- 
havior of the system on various problem parameters at 
high and intermediate T, while we support our numer- 
ical results with analytical arguments. It is well known 
that magnetic impurities get screened by the chain at 
low T forming with their neighbors an effective impu- 
rity of different spin. 24-26 Notwithstanding, as we show 
in this work, conspicuous excitations, which cannot be 
described by the effective impurity picture, survive at 
low T; this holds for a single impurity in a finite sys- 
tem or plausibly for a finite concentration of impurities 
in the thermodynamic limit. Furthermore, the screening 
of the impurity triggers Kondo-type, many body, phe- 
nomena below a crossover temperature, and the chain 
becomes perfectly (insulating)transmitting at T = for 
(anti)ferromagnetic easy axis anisotropy A. 27 ' 28 This is 
an analog of the behavior of a Luttinger liquid in the 
presence of a non-magnetic impurity where the chain is 
(cut)healed for (repulsive) attractive interactions. 29 We 
seek evidence of the cutting-healing behavior of the chain 
using the thermal conductivity as a probe, which is ideal 
for this purpose since the only energy current relaxation 



mechanism is induced by the impurity. 

From the above it is obvious that the role of the ther- 
mal conductivity in these systems is exceptional and 
twofold. Not only is the thermal conductivity a very use- 
ful theoretical tool but it is of high experimental and tech- 
nological interests as well. Thus, our results could serve 
as qualitative guidelines for forthcoming thermal conduc- 
tivity measurements in doped cuprates, e.g., Ni doping in 
the renowned SrCu0 2 and Sr 2 Cu0 3 cuprates. Further- 
more, magnetic impurities in highly heat conducting ma- 
terials could function as potential switching mechanisms 
enabling tunable heat transport and the emergence of 
numerous technological applications. Therefore, it is vi- 
tal to understand the fundamental properties of a proto- 
type model in the presence of magnetic impurities before 
trying to exploit the transport properties of these truly 
unique systems. 



II. MODEL 

In a previous work, 28 the thermal transport properties 
of a spin-| Heisenberg chain was studied in the presence 
of a spin-,!? magnetic impurity located out of the chain. 
In addition, a uniform chain with two consecutive weak 
links was discussed, which can be considered as a special 
case of a spin-i impurity embedded in the chain. In the 
present work, we deal with the generic case of a spin-S" 
magnetic impurity embedded in the chain, Fig. 1. Many 
of the conclusions and the basic ideas reached in Rcf. 28 
can be applied in the present work, as well, due to the 
single impurity character of the problem. However, sig- 
nificant differences in the physical results arise due to the 
different geometries of the problems. A brief discussion 
and comparison between the two models is given at the 
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FIG. 1: Spin-S impurity embedded in a spin-| chain. 
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end of the manuscript. 

The Hamiltonian of an AHM on a ring of L sites, where 
a spin-i operator, say the s , is substituted by another 
one of different spin, say S, is given by 

L-2 

H = J2 h ' + J ' ( s i + s i-0 ■ § . {h = M ■ si+i) , (1) 
l=i 

see Fig. 1 as well. In Eq. (1), J is the antiferromag- 
netic exchange coupling — with energy units — between 
spin- 4 operators, J' is the host-impurity coupling, while 
the tilded operators s, S denote the vector operators 
(s x , s v , As z ), (S x , S y , AS Z ), respectively. In addition, we 
work in a system of units where the lattice constant a, 
and the Planck and Boltzmann constants are a,H,kB = 1; 
yet numerical results are presented for J = 1. 

The spin f , energy f current operators are determined 
by the continuity equation d t Oi+VJi = 0. Taking 0[ to 
be the local spin, energy operators and Ji to be the local 
spin, energy current operators respectively, we arrive at 

L-2 

f = ^2 J ( Sl X S ' + ■ iz + J ' [( SL -! ~~ Sl ) X S ] ' 6z , (2) 
1=1 

where e z is the unit vector along the z-axis, and 

L-2 

f = J2SI ' ( §i +! X ^"O + JJ ' Sl ' ( §2 X ^) 

1=2 

+ JJ'SL-X • (S X § L _ 2 ) + J' 2 S • (h X S £ _i) . (3) 

Within linear response theory, the real part of the spin 
a', thermal k' conductivities are given by 30-32 

(/(w) = 2irD s 5(uj) + <j(<jS) , re'(cj) = 2ttD e 5(ljj) + k(lj) , 

where the corresponding spin D$, energy De stiffnesses 
denote the presence of ballistic transport in the system. 
The regular components cr, n of the corresponding spin, 
thermal conductivities are given by 

7T 1 — e"^ v-^ 
ct(w) = — > p n \(n\f\m)\ 5(uj mn - u), (4) 

Li UJ 

= — T /] Pn\(n\f\m)\ 2 S(uj mn - u), (5) 

where e„ are the eigenvalues and \n) are the eigenstates 
of the Hamiltonian (1), p n = exp(—(3e n )/Trexp(—f3H), 
u mn = e m - £«, and P = 1/T. 

While in the pure AHM it is clear that De is finite for 
any value of the easy axis anisotropy A, 20-22 there is an 
ongoing debate on the behavior of the spin stiffness. 33-41 
Nevertheless, a single impurity renders ballistic transport 
incoherent and both D$, De vanish. 42 Thus, the static 
component of the transport quantities is given by the 
resistive dc spin, thermal conductivities obtained by the 
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FIG. 2: (color online). The dc value of the thermal T 2 Kdc, 
spin Todc conductivity is shown as a function of the host- 
impurity coupling J'/J, obtained via ED for L — 15, S = 1, 
A = 1, and f3 -> 0. 

zero frequency limit of the regular components (4), (5), 
viz., a dc = a(u) -> 0), Hdc = k(uj -> 0). 

To numerically study transport quantities for systems 
with a Hilbert space dimension up to T> ~ 10 4 , at high 
temperatures, we employ the exact diagonalization (ED) 
technique. The 5-peaks at the excitation frequencies are 
binned in windows 8u> — 0.01 while we introduce an 
additional broadening r\ = 0.03 using the well-known 
Kramers-Kronig relations. For T> > 10 4 , we use the mi- 
crocanonical Lanczos method 43 (MCLM) at high tem- 
peratures and the finite temperature Lanczos method 44 
(FTLM) at finite temperatures; yet we keep the same 
additional broadening. Lastly, all results are obtained in 
the SfL.i = subsector. 



III. NUMERICAL RESULTS 

A. Spin-1 impurity 

First, we would like to address the S — 1 impurity case 
since this may be the most appealing magnetic impurity 
doping for experiments. In Fig. 2, we present results 
for the dc spin adc, thermal Kdc conductivities, while in 
Figs. 3, and 4 we present the dynamical thermal con- 
ductivity of the isotropic (A = 1) Heisenberg chain. A 
wide range of host-impurity couplings, J'/J = 0.4 — 2.0, 
is shown at high temperatures, j3 — > 0, for various lat- 
tice sizes, L = 15 — 21. Results for L = 15 are obtained 
via the ED technique, while results for L > 15 are ob- 
tained by employing the MCLM. In order to eliminate 
unimportant for this discussion prefactors of the trans- 
port quantities, we plot either the normalized thermal 
conductivity, R{uj) — k(oj) j kq where the normalization 
kq is given by kq — J K(ui)duj, or the non-trivial at high 
temperatures Tad Cl Ta{uS), T 2 Hd c and T 2 k(uj). Finally, 
for the discussion of the lattice size scaling we plot the 
scaled thermal conductivity k(uiL)/L. 28,42 
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FIG. 3: (color online). The frequency dependence of the nor- 
malized and scaled thermal conductivity k(ujL)/L is shown 
for L — 15,19,21 and J'/ J — 0.8; in addition a Lorentzian 
fit is depicted. Results are obtained for S = 1, A = 1, and 
ft — > 0. L = 15 results are obtained via ED and L > 15 via 
MCLM. 



FIG. 4: (color online), a) The frequency dependence of 
the normalized and scaled thermal conductivity R(ujL)/L is 
shown for L = 15, 19,21 and J' / J — 1.5. b) The frequency 
dependence of T 2 k(loL)/ L is shown for L = 15 and two host- 
impurity couplings J 1 J J — 0.4,2.0. Results are obtained for 
S = 1, A = 1, and /3 — > 0. L — 15 results are obtained via 
ED and L > 15 via MCLM. 



For a very weak or very strong J', a severe reduction of 
Kdc from its maximum value, which occurs at J' = 0.8J, 
is illustrated in Fig. 2. On the contrary, the behavior of 
(Jdc is qualitatively different from the one of Kdc- Besides 
the qualitative difference, it is rather impressive that Ode 
hardly changes in the wide range of J' shown in Fig. 
2. Although it is reasonable that a will be less sensi- 
tive to the effect of the single impurity due to the bulk 
scattering — [H, f] ^ in the pure model — its rigidity is 
still surprising. 

In Fig. 3, R(ujL)/L is depicted for J' = 0.8 J and a 
Lorentzian fit as well, signifying the Lorentzian behavior 
of k, k(u>) = Kdc/W + ( wr£ ) 2 ] with r e the scattering time. 
The thermal conductivity retains its Lorentzian shape, 
approximately, in a range of values of the host-impurity 
coupling, J'/ J ~ 0.8 ± 0.2. As previously proposed, 28 a 
Lorentzian form of k(u) is an indication of a weak per- 
turbation, while on the contrary a non-monotonic form 
implies that the system has flown to the strong perturba- 
tion regime. The Lorentzian behavior is retained, with 
an almost constant r e , for temperatures as low as the 
limiting FTLM temperature for finite size systems Tf s 
below which FTLM results are not reliable; 44 we esti- 
mate Tf s / J ~ 0.3. Although Tf s is sufficient to reveal 
the cutting of the chain for strong perturbations, 24 for 
weak perturbations a lower Tf s is required. 28 To obtain 
an impression of TV S) consider that for systems like the 
Sr 2 Cu0 3 compound with J/k B « 2500 K, T fs « 750 
K. However, defects may reduce J of doped samples 26 
bringing Tf s closer to room temperature. 

For extreme values of the coupling J', either strong 
or weak, k(lu) exhibits a strongly non-monotonic behav- 
ior, Figs. 4(a), and 4(b), indicating that the system cou- 
ples strongly to the impurity in both cases. Whereupon, 
the low frequency behavior (Luj/ J < 7r), corresponding to 
an open-like chain, 45 is similar for a weak and a strong 
J', Fig. 4(b). On the other hand, the high frequency 



behavior of k(w) for a weak and a strong J' is strik- 
ingly different due to the emergence of a conspicuous 
secondary structure at high frequencies. The frequency 
of this structure shifts with J', indicating that its ori- 
gin is local excitations of the impurity. Moreover, the 
larger the J', the more weight is accumulated at this 
structure, which becomes the prevalent contribution to 
k(u>) for fairly strong couplings, despite being only a 1/L 
effect. 



B. Spin-S* impurity 

The fact that the contribution of a 1/L effect to k(lu) 
surpasses the bulk contribution may seem quite bizarre, 
however, we can comprehend the origin of this effect from 
the analytical expression of the sum rule of k(cj). Since 
the energy stiffness De vanishes in the presence of im- 
purities the sum rule of the thermal conductivity k'(u>) 
is equal to the zeroth moment (kq) of the regular part 
k(lu). Starting from Eq. (5) and taking the high temper- 
ature limit the sum rule of the thermal conductivity can 
be written as the thermodynamic average of the square 
of the energy current operator 46 



K (oj)dcj= ni ^(ff), with, (O) 



Tre~^ H 6 
Tre-P H ' 



( 6 ) 

and 7T/3 2 (j e j e ) / L — kq. Taking the infinite temperature 
limit (/3 —¥ 0) in the evaluation of the thermodynamic 
average, Eq. (6), one arrives at 
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1 + 2A 2 (J 



J^2 
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2 J + J' 



(7) 



B 2 = \ J' 2 S{S+l) is the characteristic spin impurity de- 
pendence. The bulk contribution, ~ J 4 , to kq is equal to 
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FIG. 5: (color online). The frequency dependence of the 
thermal conductivity T 2 k(lj) is shown for two impurities 
S = 3/2,2 with J'/J = 0.6,0.5 and L = 14, 13, respectively. 
Inset: T 2 Kdc as a function of the coupling J'/J for S = 3/2, 2. 
Results are obtained via ED for A = 1, /3 — >■ 0. 




FIG. 6: (color online). The frequency dependence of the spin 
conductivity Ta(uS) is shown for various perturbations: S = 
1,3/2,2 with J'/J = 1.0 and L — 15,14,13, respectively, as 
well as the case S = 1, L = 15 with J'/J = 1.5. Results are 
obtained via ED for A = 1, /3 — > 0. 



J' 4 (we omit the - {J J') 2 
J* = JiVf ofe~, 3 -,Y - For one 

2 



the impurity contribution, 

terms), for a coupling J 1 r — „ — „ y 4 s ^ s+1 - j 

thing, for a finite system studied via ED of L = 13, S 
we have J* ~ J. Thus, k(oj) will exhibit resonant modes 
for strong perturbations, at frequencies us ~ J' — at least 
at high temperatures. 

Another important conclusion that arises from the 
evaluation of J* is that the larger the S the smaller the 
J*, since J* ~ l/\/~S. Thus, the local excitations of 
the impurity dominate even more easily over the bulk 
contribution for larger S. As a matter of fact, we did 
not obtain a Lorentzian k(lj) for any S > 1 by adjust- 
ing the host-impurity coupling, implying that any S > 1 
constitutes a strong perturbation for k(uj) while 5 = 1 
seems to be quite unique. As far as o{u>) is concerned, 
its high frequency regime is dominated by the impurity 
local excitations, similarly to k(ui), but its low frequency 
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FIG. 7: The energy spectrum of the isotropic 3-spin system, 
Eq. (8). 



a(u) — > 0) regime remains virtually unaffected for higher 
S. This is a striking difference between and k^c where 
the former shows a surprising rigidity to the influence of 
the impurity (large S, strong/weak J') while the latter is 
severely reduced in the strong perturbation regime. The 
above arguments are summarized in Figs. 5, and 6 where 
the frequency dependence of the thermal, spin conductiv- 
ity is shown respectively. In Fig. 5, T 2 n{uS) is shown for 

5 = 3/2,2 and J'/J = 0.6,0.5, namely, the correspond- 
ing couplings for which is maximum, Fig. 5 inset. For 
the case of the spin transport we present Ta(uj) in Fig. 

6 for various moderate-strong perturbations. 



C. Lattice size scaling 

Let us now turn our attention to the lattice size scal- 
ing. A Lorentzian k(lu) trivially obeys a universal L scal- 
ing, with the size independent quantity to be k(u)L)/L, 
since Kdc, T e ~ L, Fig. 3. 28 On the other hand, for strong 
J', k(lu) exhibits two L-scalings; the prominent impurity 
contribution at u) ~ J', which is 0(1), does not scale 
with L, Fig. 4(a), while the low frequency bulk contri- 
bution obeys the proposed scaling for any J'. Moreover, 
in accord with the evaluation of J* (Eq. (7)), Fig. 4(a) 
shows that the contribution of a single impurity dwindles 
with respect to the bulk contribution with increasing L, 
becoming negligible in the limit L — > oo. However, con- 
sidering the thermodynamic limit and a finite but dilute 
impurity concentration cj — which actually would be a 
more pragmatic approach to a real system — one could 
plausibly assume that the high frequency structure will 
be present, and similar scaling behaviors would hold with 
the substitution 1/L —¥ cj. 



IV. STRONG COUPLING LIMIT 

Interesting conclusions can be reached in the strong 
host-impurity limit, J' — ¥ oo, and are rather useful to 
understand the high frequency behavior of k(uj), a(u>). 
Starting from Eq. (1) for the isotropic point (A = 1) and 
taking J = 0, we end up in the Hamiltonian of a 3-spin 
system 



H = J'( Si _i + si)-S 



(8) 
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FIG. 8: (color online), a, b) The frequency dependence of the 
thermal conductivity Tk(lu), spin conductivity <t(oj) respec- 
tively for an S = 2 impurity with J' = 2J, L = 19, A = 1 
and two temperatures T/J = 0.7, 2. Results are obtained via 
FTLM. 



FIG. 9: (color online). The frequency dependence of the nor- 
malized thermal conductivity k(lj) for S = 1, J' = 2 J, L = 19 
and A = ±0.5 at T/J — 0.7,50. Results are obtained via 
FTLM. 



Exploiting the rotational symmetry and the limited de- 
grees of freedom of H, we perform an analytical diag- 
onalization into different 5 z -total subsectors. The low- 
est/highest \S Z \ = 5 + 1 subsectors are T> = 1 Hilbert 
spaces. The second lowest/highest \S Z \ — 5 subsectors 
are V = 3 Hilbert spaces. The rest 25 + 1 5 Z subsectors, 
with 5 = 5 — 1, arc T> = 4 Hilbert spaces. The energy 
spectrum of the 3-spin system is shown in Fig. 7. 

Considering the local energy current operator f — 
J /2 S • (§i x §£_i), we can see that its matrix elements 
vanish between non-zero eigenvalues, and, consequently, 
only transitions between zero and non-zero eigenvalues 
survive; this holds apart from the isotropic point as well. 
Particularly, for the isotropic Heisenberg model, the only 
non- vanishing transitions correspond to an energy differ- 
ence Se = +</'. Thus, for A = 1, k(lu) will exhibit only 
one sharp peak at high frequencies located at oj ~ J' 
and it will be independent of 5. Transitions correspond- 
ing to Se — ±J' are between elevated eigenstates, Fig. 
7, thus the peak at u> ~ J' is expected to diminish with 
decreasing temperature and finally to vanish as the sys- 
tem flows towards its ground state. These selection rules 
do not hold for the spin transport where there are more 
allowed transitions. These transitions involve the ground 
state e g as well and consequently resonant peaks will be 
present at low T. Lastly, for A^l there are transitions 
of f which do not necessarily correspond to Se = ±J' 
and yield a more complicated high frequency structure 
for k(u>). 

The previous analytical arguments are verified in Fig. 
8 where results are obtained via FTLM for L = 19 and 
various perturbations. The frequency dependence of the 
thermal conductivity Tk(uj) Fig. 8(a) and the spin con- 
ductivity cr(w) Fig. 8(b) is shown, for 5 = 2, J' = 2J, 
A = 1 at T/J = 0.7,2. 

Starting with Fig. 8(a) we can see that at high tem- 
peratures the prevalent contribution to k(u>) comes from 
the operator f yielding a prominent peak at uj = J' 



independent of 5 (compare Figs. 4(b) and 8(a)). As the 
temperature decreases and the system flows to its ground 
state the peak at uj = J' decreases gradually and eventu- 
ally vanishes. Note that the high frequency leap in Fig. 
8(a), present at T/J = 0.7, emerges from cx J J' terms 
of the f current, Eq. (3). In contrast to the thermal 
transport, the local excitations of the spin current in the 
reduced system, Eq. (8), involve the ground state giving 
a sharp peak at w = J' (5 + 1) which does not vanish 
with decreasing temperature, Fig. 8(b). 

In Fig. 9 we present the normalized thermal conductiv- 
ity k{uj) for 5 = 1, J' = 2 J, A = ±0.5 at T/J = 0.7, 50; 
results are obtained via FTLM for L = 19. First, in con- 
nection with the previous arguments, we can say that the 
complicated high frequency structure of k(uj) for A ^ 1 
indicates that the single excitation at uj = J' is a unique 
property of the isotropic Heisenberg model. Second, wc 
focus on the low frequency part of k(lu) and particu- 
larly on the stark difference for A ^ at low temper- 
atures. The chain exhibits cutting(hcaling) behavior for 
A > 0(A < 0) for an 5 = 1 impurity embedded in the 
chain which was previously reported for an 5 = J impu- 
rity out of the chain. 27 ' 28 At high temperatures the curves 
for A = ±0.5 are one on top of the other. As the tem- 
perature decreases, k(ui) for A — —0.5 tends to obtain a 
more Lorentzian-like form, characteristic of the weak per- 
turbation regime. On the contrary, for A = ±0.5, k(u>) 
obtains a strongly non-monotonic behavior with decreas- 
ing temperature resembling the thermal conductivity of 
an open chain and signifying the flow to the strong per- 
turbation regime. Note that cr(ui) exhibits a similar low 
frequency behavior for A = ±0.5 which could be plau- 
sibly read as evidence of a finite D$, yielding a ~ 6(ui) 
contribution to u'(oj), in the pure AHM for |A| < 1. 

It is also interesting to present the temperature de- 
pendence of the dc value of the thermal conductivity k<j c 
itself since this is a directly measurable quantity in ex- 
periments. In Fig. 10, TKdc(T) is shown for four charac- 
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FIG. 10: (color online). The temperature dependence of 
TK, dc {T) for 5 = 1, J'/ J = 0.8,2 and S = 2, J'/ J = 1,2; 
results are obtained for A = 1, L = 19 via FTLM. 



teristic perturbations and A = 1. The main conclusion 
is that the maximum of Kdc occurs at different temper- 
atures for different perturbations. Although for 5 = 1 
the host-impurity coupling quantitatively only affects the 
behavior of Kdc this is not the case for other impurities 
like the 5 = 2 impurity. Generally the maximum of Kdc 
may occur at different temperatures even for the same 
impurity if this is coupled to the chain with different 
host-impurity couplings J'. 

Finally, we would like to address the screening of the 
impurity by the chain which has been discussed previ- 
ously in the literature. 24-26 The ground state of the 3-spin 
system e g is (25 + l)-fold degenerate, Fig. 7, implying 
that the impurity with its neighbors form an effective spin 
5 at low energies. Thus, for a strong J' one can assume 
that the degrees of freedom of the system at low energies 
will be described by states of the form \e g )®\4>L-3) , where 
the pseudo spin 5 couples with the rest of the chain with 
an effective coupling J. J can be evaluated considering 
the matrix elements of the operator hi at low energies. 47 
For an 5 = 2 impurity we obtain the effective coupling 
to be ferromagnetic, J = — 0.25J, while the larger the 5 
the weaker the J. Thus, for a single impurity in a finite 
system, or a finite concentration of impurities in the ther- 
modynamic limit, the picture of the effective spin clearly 
fails to describe the whole frequency range since the weak 
J, | J| <C | J' |, cannot reproduce the high frequency, con- 
spicuous, excitations yielded by a strong J' at uj ~ J'. 



V. DISCUSSION AND CONCLUSIONS 

It is worthwhile to devote a few lines to contrast the 
basic points of the present model, where the impurity is 
embedded in the chain (IEC), with those of the model 
studied previously with the impurity located out of the 
chain (IOC). 28 Stark differences arise in the transport 
properties of the two models due to the position of the 
impurity and the way it couples to the pure system. For 
instance, a weak host-impurity coupling J' is a strong 



perturbation for the IEC model, Fig. 4(b), while taking 
J' = in the IOC model we end up in the pure AHM. 
Note also that it would not be accurate to perceive the 
difference \J — J'\ as a perturbative parameter, simply, 
because the case J' = J docs not correspond to the min- 
imal perturbation of the heat transport. In addition, for 
different impurities the maximum Kd c occurs at different 
J' which does not necessarily correspond to the maxi- 
mum of o~dc, Figs. 2, 5 inset. For spin-5 impurities with 
5 > 1 k(oj) exhibits a non-monotonic form. The absence 
of a Lorentzian k(uj) for 5 > 1 and any host-impurity 
coupling J' is an indication that 5 > 1 impurities con- 
stitute a strong perturbation for the heat transport of 
the Heisenberg model. This is in sharp contrast to the 
behavior of k(uj) in the IOC model which obeys a uni- 
versal scaling with the B" parameter, at least for weak- 
intermediate J'. 

Another significant difference between the two models 
arises from the absence of a oc J' 2 term in the energy cur- 
rent of the IOC model. As we have shown in this work, 
the j e term accounts for the prominent high frequency 
excitations which become the prevalent contribution to 
k(uj) for strong perturbations. Similarly, the high fre- 
quency behavior of o~(uj) for strong perturbations is dif- 
ferent in the IEC and the IOC models due to the absence 
of a perturbative spin current term oc J' in the latter. 
From all the above one could conclude that the magnetic 
impurity in the IEC model is a much stronger perturba- 
tion for a Heisenberg chain than in the IOC model. 

To summarize, studying the thermal k and spin a 
conductivities of the ID, spin-i, AHM with an embed- 
ded spin-5 impurity at finite temperatures T, we have 
reached the following conclusions: (i) An 5 = 1 impurity 
can be considered as a relatively weak perturbation for 
some host-impurity couplings J', in contrast to 5 > 1 
impurities which have a more drastic effect on k. Fur- 
thermore, the difference between Ode and Kd c is remark- 
able since the former shows an impressive rigidity under 
the influence of the impurity, (ii) k(oj) obeys a univer- 
sal scaling k(uL)/L for weak and intermediate J', while 
a strong J' triggers the emergence of impurity local ex- 
citations ruining the L-scaling of k(oj). (iii) We have 
demonstrated the origin of these resonant modes for a 
large J' and a finite system — their position, their mag- 
nitude, and their temperature behavior as well, Eq. (7), 
Fig. 8. (iv) The presence of these sharp excitations at 
low T makes the picture of the effective spin 5 = 5 — 1 in- 
sufficient to describe the finite frequency transport prop- 
erties, Fig. 8 — at least for a single impurity in a finite 
system or plausibly for a finite concentration of impuri- 
ties in the thermodynamic limit, (v) Finally, we observe 
the cutting-healing behavior of the chain according to the 
sign of A as it was previously reported for magnetic im- 
purities out of the chain, Fig. 9. 27,28 
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